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Abstract 

We examine the phenomenology of the CP-odd scalar ^4° of two-Higgs-doublet 
models. We explore the parameter space determined by triviality bounds and iden- 
tify the regions where the ^4° can be detected at the LHC in each of the following 
modes: the inclusive two-photon decay mode, the I 77X mode from tiA° produc- 
tion and the A — > Zh channel with the subsequent decay of the CP-even scalar h 
to two photons. We find that, while the l^j^/X mode is of limited usefulness, the 
other two modes can give viable signals in fairly large, and complementary, regions 
of parameter space. 
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1 Introduction 



Despite the remarkable experimental success of the Standard Model of electroweak inter- 
actions, there is no evidence, so far, in favor of the minimal one-doublet scalar sector of 
the model. It is therefore essential that alternatives to it be considered. The simplest 
extension of the scalar sector of the Standard Model that naturally embodies certain in- 
dispensable features, such as the smallness of p — 1 and the suppression of flavor-changing 
neutral currents, while admitting the possibility of new observable phenomena, is the 
(CP-conserving) two-Higgs-doublet model. 

In the Standard one-doublet Model, the Higgs mass is constrained to lie between 57 
and 800 GeV approximately. The lower limit is set by direct searches at LEP []]], while 
the upper bound comes from theoretical considerations such as the triviality of theories 
with fundamental scalars 0,0,0. The numerous phenomenological investigations of the 
one-doublet Standard Model indicate that it may be possible to explore this entire 
range at future colliders such as LEP-II and the LHC So we may ask about the potential 
of these colliders to detect alternative models. For a general two-Higgs-doublet model, 
a first step towards answering this question was made in ||, where the bounds on the 
parameters of the model that result from triviality considerations were derived. These 
bounds were obtained in the context of perturbation theory and so their validity is of 
a rather qualitative nature. To be definite, we will focus our interest on the region of 
parameter space defined by these triviality bounds. However, the general discussion we 
offer can give a fair idea of how to extend our results beyond this region. 

In the analysis of ref. ||, following the spirit of |§, the term 'triviality' was employed 
in the following sense: The running scalar self-couplings of the Two-Doublet Model are 
expected, at least in perturbation theory, to develop a Landau pole at a finite momentum 
scale. Consequently, the theory can only stand as an effective low-energy theory valid up 
to some finite cutoff A, beyond which new phenomena emerge. The calculation of physical 
quantities in the effective theory will thus be accurate up to terms of order pf/A 2 , M 2 /A 2 , 
where pi are typical momenta of the process under consideration and Mj are the masses 
of the particles of the theory. In ref. |J, the effective theory was defined to be valid if all 
masses satisfied 

A - 2tt v ; 

Thus, given a set of parameters of the model, that is, masses and couplings, a cutoff was 
defined by 

A = 2vr max M,- (2) 

3 

and the following consistency requirements were made: (i) No coupling should develop a 
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Landau pole at a scale less than A, and (ii) the effective potential should be stable for all 
field values less than A. From these two conditions follow the 'triviality bounds' reported 
in H and used in the present study. 

It is thus our intention to find out how much of the parameter space allowed by 
triviality can be explored at present and future colliders. In this paper we concentrate 
on the phenomenology of the CP-odd scalar (which we denote A ) of the Two-Doublet 
Model at the LHC, and identify those regions of parameter space where direct observation 
of this particle is possible. 

We examine three of the most promising signals for the A . First, we discuss the 
inclusive two-photon channel. As is the case for the Standard Higgs, a light A decays 
predominantly into bb pairs and thus its detection at a hadron collider has to rely on rare 
modes, such as the two-photon mode. This signal usually consists of a large number of 
events and, furthermore, it may be useful even beyond the so-called 'intermediate mass' 
range, because the A does not couple to a pair of weak gauge bosons at tree level. On 
the other hand, it is well known that there are very large backgrounds to this signal, 
whose successful containment places severe requirements on the design of the detector 
0. Consequently, even a large event rate may not constitute a clean signal. For this 
reason, we also examine the process in which the A is produced in association with a 
tt pair and decays into two photons while one of the t/t is tagged by its leptonic decay 
f|||. The lepton tag ensures a much cleaner signal, but the event rate is so small as to 
render this mode only marginally useful. Finally, we turn to the discussion of the decay 
A — > Zh, where h is a CP-even scalar, with the subsequent decay of h in two photons, 
and demonstrate that it gives a very clear signal in a region of parameter space where 
none of the two other modes considered can be useful. To our best knowledge, this process 
has not been considered before in this context. Our results indicate that it may provide 
an excellent way of directly observing the CP-odd scalar A . 

Phenomenological studies of two-Higgs-doublet models have been carried out mostly 
in the context of the Minimal Supersymmetric Model (MSSM) |10|,|ll],|12|,|13|,|l|,|15[ ; if the 
supersymmetric partners are heavy enough, then the low-energy spectrum and the formal 
structure of the couplings of this model are exactly the same as in the more general 
non-supersymmetric case. Supersymmetry, however, imposes relations among the various 
masses and couplings, which reduce the number of independent parameters and allow for 
more definite predictions. Moreover, the parameter space usually investigated in studies 
of the MSSM is motivated by supersymmetric unification, while in the present study it 
is defined by the above-mentioned triviality bounds and, as a result, it is significantly 
different. It might be thought that the MSSM should really be contained in the more 
general two-doublet model as a special case. The models studied here, however, employ an 
exact discrete symmetry (see Section 2 and ref. flBfl ) intended to suppress flavor-changing 
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neutral currents at tree level; if this symmetry is allowed to be softly broken, as is the 
case in the MSSM, most triviality bounds no longer hold ||. In this case the parameter 
space of the MSSM is indeed a subset of that of the more general two-doublet model. 

The larger number and different range of the parameters of the Two-Doublet Model 
may result in a markedly different phenomenology from that of the MSSM. A striking 
example is the A — > Zh — > Z77 decay sequence, described in detail in Section 7. We 
will further comment on the similarities and differences between the supersymmetric and 
non-supersymmetric two-doublet models in the course of the discussion of our results. 

In the next section we review the Two-Doublet Model, emphasizing the features most 
relevant to phenomenology. In Section 3 we describe the production mechanisms of the 
CP-odd Higgs A at the LHC In Section 4 we discuss the main decay channels and 
present branching ratios of the A in order to motivate our choice of the particular signals 
we subsequently examine. Section 5 concerns the inclusive two-photon channel, while 
in Section 6 we study the signal from ttA° production. Section 7 is devoted to 

the discussion of the decay A — ► Zh. The results we present in Sections 5-7 consist 
in event-rate and significance contours in appropriate sections of parameter space, for 
integrated luminosities of 10 fb _1 and 100 fb _1 , which correspond to a low- and a high- 
luminosity option at the LHC. We have used a center-of-mass energy of ^/s = 16 TeV in 
all calculations. Finally, Section 8 contains our conclusions. 



2 The two-doublet model 

The scalar sector contains two electroweak doublets $1, $2, both of hypercharge Y — 1. 
In order to eliminate flavor-changing neutral currents at tree level, one has to impose a 
discrete symmetry |16| . There are many ways of doing this, the two most often discussed 



being the following^ JT7 , 18 



• Model I : $ x -> -$1 

(3) 

• Model II : $1 —> -$1 ; d Ri -> -dm, e Ri -> -e Ri 

{dm (i = 1,2,3) are the right-handed negatively charged quarks and em the right-handed 
charged leptons.) The Lagrangian is 

C = Ckin + £>Y — V 



: The notation here has been changed from that of ref. Q to facilitate comparison with most studies of 
the MSSM. Contact with || can be made by the following substitutions: <t>i «-> $2, Model I Model II. 
In terms of the angles a and /3 introduced later in this Section, this has the effect a —a, [3 «-> 7r/2 — (3. 
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where Ckin contains all the covariant derivative terms, V is the scalar potential and Ly 
contains the fermion-scalar interactions. The form of the latter is the following: 

• Model I 

£y = g^ku^Rj + Q^ku^dRj + g^'lu^eBj + h.c. (4) 

• Model II 

£y = 9ijlu^2 u R3 + 9ij ] QLi^i d Rj + 9ij~ l Li^ie Rj + h.c. (5) 

where qn are the left-handed quark doublets, um are the up-type right-handed quarks 
and In are the left-handed lepton doublets. The matrix g"' is understood to be diagonal. 

Thus in Model I only $ 2 couples to fermions, while in Model II $x gives mass to 
down-type quarks and leptons and $ 2 to up-type quarks. Because of the different fermion 
couplings, the two models have different phenomenology. We will be presenting results 
for each of the two models^. 

The discrete symmetry ([3D also limits the number of scalar self-couplings. The scalar 
potential is 

V = ^$I$i + ^$^ 2 + A 1 ($ t 1 $ 1 ) 2 + A 2 ($|$2) 2 + A 3 ($ t 1 $ 1 )($|$2) 

+A 4 ($I$ 2 )($£$i) + \K [{^2? + ($^i) 2 ] (6) 

Provided certain inequalities among the scalar self-couplings are satisfied |Tj| , the doublets 
$i, $ 2 acquire vacuum expectation values of the following form: 

where v% , v 2 are real and 

v\ + v 2 2 =v 2 = (246 GeV) 2 . (8) 

Note that CP is a good symmetry of the scalar sector and so we may assign definite CP 
quantum numbers to all neutral scalar states. At higher loop level, interactions with the 
fermions will invalidate this statement and mixings between CP-odd and CP-even states 
will appear due to the CP- violating phases of the Kobayashi-Maskawa matrix. These 
mixings, however, are small, and we will therefore neglect them in what follows. 

2 It should be noted that the MSSM is of type II only. 
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Table 1: Weak gauge boson and fermion couplings of the neutral scalars of the two-doublet 
model relative to the corresponding Standard Higgs couplings. 



In addition to the Goldstone bosons which become the longitudinal components of 
the W's and the Z via the Higgs mechanism, the spectrum of the scalar sector con- 
tains two neutral CP-even scalars, denoted by h, H, one neutral CP-odd scalar, A , 
sometimes referred to as a pseudoscalar due to its 7 5 coupling to fermions, and two 
charged complex conjugate states H ± . It is customary to introduce two mixing angles: 
(3 (= arctan(t> 2 /f i)) rotates the CP-odd and the charged scalars into their mass eigen- 
states while a (—7i/2<a< 7r/2) rotates the neutral scalars into their mass eigenstates. 
There are six independent parameters in the scalar potential (after fixing v through (|8|)) 
which can be taken to be a, (3, Mh±, Mh, Mh, M^o. (In contrast, in the MSSM there are 



only two.) To these one should add the top quark mass, expected to lie in the range [2(| 



100 GeV £ M t £ 200 GeV (9) 

but otherwise unspecified. 

In Table 1 we display the couplings of the neutral scalars to vector bosons and fermions 
for Models I and II. The values shown are relative to the Standard Model Higgs couplings; 
V denotes collectively the weak gauge bosons, u the up-type quarks and d the down-type 
quarks and leptons. We also display below the couplings of A to a scalar and a vector 
boson, since they are important in the determination of the partial widths of A . 

1 



Ca°sv = 2 VV 2 + 9' 2 [A & 1 (H sin(/3 -a)-h cos(/3 -a))]--g W+(A° &> H~) + h.c. 

(10) 

where g , g are the U(l)y and SU(2) gauge couplings respectively. The full set of Feynman 
rules can be found in ref. flT7| . 
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Of crucial importance to the phenomenology of the CP-odd Higgs A is the absence of 
a tree level coupling between it and a pair of weak gauge bosons. This deprives us of a very 
clean signal as well as of an important production mechanism for masses of a few hundred 
GeV. It is also evident from Table 1 that the couplings of the Higgs bosons to fermion- 
antifermion pairs can be significantly enhanced or suppressed over those of the Standard 
Higgs depending on the values of (3 and a and the choice of Model. Similar suppressions 
can occur in the couplings of the A to the other neutral scalars (see eq. (0)). Unlike 
the situation in the MSSM, these couplings are independent of the scalar masses. This 
freedom has important consequences for phenomenology, as will be demonstrated in later 
sections. 

The strongest experimental constraints on the two-doublet model arise from the recent 
CLEO bound on BR{b — > s-y) []2TJ] and from LEP data on the Z — > bb decay width 



22| . These measurements impose a combined bound on (3 and M#± [p3|,P4| . There 



are uncertainties p3) about the precise value of these bounds, and we will not include 



them in our analysis. From refs. [p3p4| , however, it is clear that light charged scalars 
(M H ± ^ 50 GeV) and low (3 {(3 & 15 — 20°) are strongly disfavoured, with the restrictions 
becoming tighter for larger top quark masses. 



3 Production mechanisms 

As was emphasized earlier, the coupling of the A to intermediate vector bosons is induced 
only at loop level and is therefore small. So, in contrast to the case of the Standard 
Higgs, production of the A via vector boson fusion or in association with a W or a Z is 
insignificant. Thus, at a hadron collider, we have to rely on the following mechanisms: 

(a) Gluon fusion gg — ► A 

(b) Associated bb~A° production qq, gg — > bb~A° 

(c) Associated it A production qq, gg — > it A 

The cross sections for processes (a)-(c) are identical to the corresponding rates in the 
MSSM, since they only depend on M^Mj and tan/3. We discuss them here for com- 
pleteness and for future reference. Our results are in agreement with ref. |T3| . 



3.1 Gluon fusion 



This is usually the dominant mode. The reaction proceeds via a quark loop |25]. The 



contribution of the top quark is generally the most significant, but, in Model II and for 
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large (3, the bottom quark contribution can become equally important or even dominate, 
especially for small M A o. We have taken into account leading QCD corrections as follows: 
in ref. |26| the QCD corrections to the production of a pseudoscalar have been calculated 
in the limit of a very heavy quark flowing in the loop. The major part of these corrections 
can be cast into the form 



a = cr 



1 + (6 + 7T 



2 a s (M A o) 



7T 



:iv, 



in analogy to the case of a CP-even Higgs |10|- cr and a are the lowest-order and QCD- 
corrected cross-sections respectively. For our purpose, this is a valid approximation only 
if M40 < 2M t and only insofar as the bottom loop contribution, to which this correction 
does not apply, is negligible compared to that of the top quark. With hindsight, we may 
say that the second qualification is not worrisome: in the regions where the bottom quark 
contribution is significant, namely large (3, small M A o, the rates for the signals we will be 
interested in are small anyway. This may not be so for /3 very close to 90°, but triviality 
bounds almost always exclude such values of (3. The QCD corrections to arbitrary mass 
pseudoscalar production from gluon fusion have not been computed. Consequently, we 
shall only employ (01]) if M A o < 2M t . 

Figures la, lb show the production cross-section of the CP-odd Higgs through gluon 
fusion as a function of its mass for Models I and II respectively. The effect of the bottom 
quark loop for large (3 in Model II is evident. The top quark mass was taken to be 
150 GeV. The curves for different M t are similar; they always peak at M A o = 2M t , as 
long as j3 does not take on values extremely close to 90° in Model II. The value at the 
peak is lower for higher M t . Here and throughout this study we employed the EHLQ 



structure functions, Set 2, with Ajjg = 290 MeV p7 . 



3.2 Associated bbA° production 

This mechanism can be important for large (3 in Model II when the bbA° coupling becomes 
strong. For M A o 3> M b , this process can be approximated by bb fusion using the 
bottom quark distribution functions within the proton. The bb fusion cross-section is 
proportional to the partial width T(A° — > bb). QCD corrections to the latter have been 
included in our calculation as we explain in Section 3. No other QCD corrections to this 
process have been estimated. 

In Figures 2a, 2b we show the production cross-section of A from bb fusion as a 
function of the A mass for Models I and II. It can be seen that, in the case of Model I, 
this cross-section never exceeds the gluon fusion cross-section. In contrast, this mechanism 
dominates in Model II for very large values of (3. 
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3.3 Associated ttA° production 

This process always gives a smaller rate than gluon fusion. It is interesting, however, in 
its own right, because it is possible to obtain a relatively clear signal if one can tag on a 
lepton from the decay of the top or antitop. This will be further discussed in Section 6. 



We have computed this cross-section using the spinor techniques of ref. [29|. QCD 
corrections have not appeared in the literature and are not included. In Figures 3a, 3b 
we display the cross-section for the process pp — > it A as a function of the pseudoscalar 
mass and for the set of values of (3 used in the previous graphs. The top quark mass is 
120 GeV and 180 GeV respectively. There is only a slight drop in rate with increasing 
top quark mass, as the smaller parton luminosities at high mass are compensated by the 
enhanced Yukawa coupling. There is no distinction between Models I and II in this case 
because the A°ti coupling is the same in both Models (see Table 1). We also show the 
corresponding curves for the Standard Model Higgs. In all calculations we have included 
the contributions of both the gg and qq initial states. The latter is negligible for the CP- 
odd scalar but amounts to 20-25% of the gluon-initiated reaction for a Standard Higgs 
lighter than about 200 GeV (see also [P0f| ). We notice that the curves for the CP-odd Higgs 
are much flatter than those for the Standard Higgs. For (3 = 45°, the it A production rate 
for A masses less than about 150 GeV is significantly lower than the corresponding rate 
for the Standard Higgs. As we will see in Section 6, this has unpleasant phenomeno logical 
implications. 

In Figures 4a, 4b we show the total production cross-section for the A , the sum, that 
is, of the three processes discussed above. Results are shown for the two Models and for 
the same values of (3 used in previous figures. A top quark mass of 150 GeV was assumed. 
We again display the corresponding cross-section for the Standard Higgs. (Gauge boson 
fusion has not been included in the calculation of the latter.) The conclusions to be drawn 
from these graphs are the following: (i) except for very large values of f3 in Model II, the 
production cross-section increases with decreasing (3; (ii) the two Models diverge only 
at large j3; (iii) for (3 ^ 50° the cross-section is larger than that of the Standard Higgs; 
(iv) the cross-section is generally peaked at M40 = 2M t . Variation of M t shifts the peak 
accordingly and lowers the value at the peak somewhat. The other parameters of the 
two-doublet model do not affect the production rate of the CP-odd scalar. 



4 Decays and Branching Ratios 

In this Section we review the major decay modes and the corresponding branching ratios 
of the CP-odd Higgs in order to identify the signals that can be of phenomenological 
interest. In the absence of a tree level coupling to the W's and the Z, the A decays 
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mainly to fermion pairs or to states containing other scalars. Loop-induced decays can 
also be large, like the two-gluon mode, or phenomenologically interesting, like the two- 
photon decayQ. The main channels are the following: 

A -> 66, cc, tt, t + t-, W ± H T , Zh, ZH, gg, 77, Z 1 . (12) 



All decay widths have been calculated to lowest non-trivial order. Analytic formulae 
for the last three, loop-induced, decays can be found in the literature piyiT|]l3] ■ QCD 
corrections have been included where they are most significant, namely: 

(i) the qq decays (q = 6, c): Leading QCD effects have been accounted for by using, in 
the expression for the Yukawa coupling, the running quark mass at the scale of the mass 
of the CP-odd Higgs f3§; 



(ii) the 77 and Z7 decays through quark loops: Corrections have been calculated in 
ref. |33| for the 77 decay. The leading contribution can be incorporated by using the 
running quark mass at scale M40. The remaining QCD corrections are mostly within 
10% and have not been included. We have also used running masses in the computation 
of the decay to Z7; 

(iii) the gg decay: Corrections to this decay have only been computed for the case of 
a CP-even Higgs f34|. In analogy to that calculation, we used running quark masses in 
the loops, but did not include any further corrections. 

The branching ratios of the CP-odd scalar are complicated functions of the seven free 
parameters of the model and hence impossible to describe completely. We will illustrate 
some of the most prominent features by means of a few examples. 

Fig. 5 shows the branching ratios of the A as a function of its mass for a case where 
the other scalars, h, H and H^, are heavier than the top quark. The set of values of the 
parameters chosen (referred to as Set A) is the following: 



Set A : M H = 400 GeV, M h = 260 GeV, M H ± = 350 GeV, 

M t = 120 GeV, a = 30°. ^ ' 

Branching ratios are shown for two values of /3 (30° and 60°) and for Models I and II. Of 
the decay channels listed in ([12]) we do not show here the branching ratios to cc, gg, ZH 
and Zj. The first two are phenomenologically uninteresting at hadron colliders. The ZH 
channel has a behaviour qualitatively similar to that of the Zh mode, while the branching 
ratio to Z7, a channel of potential phenomenological importance, is generally very small. 
We point out the following features: (i) For M40 < 2M t , the dominant decay mode is 66, 

3 See however, for an investigation of the mode A — > ZZ in the MSSM. 
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unsuitable for detection of a Higgs at a hadron collider .0 Thus, the 77 or t + t~ decays have 
to be employed, (ii) Whereas the branching ratio to two photons in Model I is insensitive 
to j3, in Model II it is a decreasing function of (3, at least for moderate (3 values. This is 
easy to understand; the two-photon width is governed by the A°tt coupling, because of 
the top loop that mediates the decay, and so goes like cot 2 /? (see Table 1), whilst the full 
width is essentially saturated by the bb mode (~ tan 2 /3). As a result, the 77 branching 
ratio has an approximate cot 4 (3 dependence. In Model I, however, the A couplings to bb 
and it are both proportional to cot /3 and the /^-dependence drops out. (iii) The channels 
that contain scalars (Zh, W ± H T ) become dominant, or at least compete with the it 
mode, as soon as they open, (iv) Both the t + t~ and the 77 branching ratios fall abruptly 
once the it channel becomes available. In the region between the it and the scalar particle 
thresholds, detection of the A seems dubious. 

With Fig. 6 we wish to illustrate the dependence of the branching ratios on the angle 
a instead. In what we call Set B in this figure, the parameters have the following values: 



Set B : M H = 400 GeV, M h = 100 GeV, M H ± = 350 GeV, 
M t = 180 GeV, a = 30°, 



(14) 



whereas in 'Set C they are 



Set C : M H = 400 GeV, M h = 100 GeV, M H ± = 350 GeV, 
M t = 180 GeV, a = -45°. 



(15) 



(3 is fixed at 40° in all graphs of this figure. This is an example of the Zh channel opening 
up before the it does. The main points we want to emphasize here are the following: (i) 77 
and t + t~ are again the only possibilities for detection of a light A ; (ii) the a-dependence 
comes from the A°Zh coupling (see eq. (|T0D) which is proportional to cos(/3 — a). In Set B 
this coupling is not suppressed and the Zh mode is dominant even after the it channel 
becomes available. In contrast, in Set C, the Zh branching ratio is sizeable only before the 
it threshold. What is also to be noticed in this case is that the opening of the Zh channel 
does not cause any significant drop to the value of the two-photon branching ratio. As a 
result, the 77 decay remains a promising mode for A masses as large as 2M t . It should 
be noted that no such suppression is present in the A°W ± H Zf coupling (see eq. (fT0|)). 



Consequently, the two-photon branching ratio necessarily drops once the W^H^ channel 
is open, independently of the value of a. 

Although the particular shapes of these curves change for different values of the Two- 
Doublet Model parameters, the general features described above persist. We are thus led 



4 (3 = 30° in Model II illustrates the case of the gg branching ratio (not shown) becoming dominant 
just below the top threshold. 
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to examine the observability of the CP-odd Higgs in the two-photon mode (as the most 
promising alternative for a light A , but potentially useful even for an A of mass as large 
as 2M t ) and in the Zh mode (which, unless it is suppressed, covers a range of higher 
masses Mao). We have not further considered the t + t~ modeQ. 



5 The inclusive two-photon channel 

This mode has been extensively studied in the literature for the CP-odd scalar of the 
MSSM [|H],|rj.[rj|.[nj . Since the parameter space under investigation here is qualitatively 
rather different, as explained in the Introduction, an independent study is necessary. 
We shall briefly compare our findings to those for the MSSM, to the extent that this is 
possible, at the end of this Section. 

To assess the observability of a Higgs signal in this mode a reliable calculation of 
the reducible and irreducible backgrounds is required. This, in turn, presupposes a good 
understanding of the detector parameters and a full Monte Carlo simulation. We did 
not make such an effort here, but only tried to approximate a realistic situation as much 
as possible, within the framework of simple parton-level calculations. Our procedure of 
estimating the backgrounds, which are the same as in the Standard Model, is sketched 
below. 

At lowest order, the direct two-photon background arises from the Born process 

qq -> 77 (16) 

The gluon fusion process 

99 -> 77, (17) 



through the 'box' diagram, gives a contribution of the same order of magnitude as (|16|) . 
We calculated the cross-sections for processes ([IB]) and fll7D using the following kinematical 
cuts to simulate the detector acceptance: 

|rf | < 2.5 ; pi > 20 GeV (18) 

In the calculation of the 'box' diagram we used five light flavours and a top quark of 
150 GeV. Variation of the top quark mass in the range 100 GeV<M t <200 GeV produces 



5 See jn| for a discussion of this channel in the MSSM. The r + r branching ratio for a light A is 
of the order of 10% in both the MSSM and the Two-Doublet Model; consequently, the findings of [(To) 
should apply here too. 
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negligible change in the cross-section. A next-to-leading order calculation of the di-photon 



background p5| , |36| shows that bremsstrahlung processes, where one or more photons are 
radiated off a jet, largely dominate over the lowest order ones. However, suitable isolation 
cuts can reduce this higher-order background to about half the sum of the Born and 



'box' contributions. We therefore multiplied the calculated cross-section of processes fllq) 
and fllTD by a 'K-f actor' of 1.5 in order to 'convert' to the next-to-leading-order result. 

The background from events where jets are faking photons can be significant. We 
assumed (following [ |37| . |3"8|j ) that it will be possible to achieve such a jet rejection efficiency 



as to reduce these backgrounds to well below the prompt di-photon level. In fact, even 
so, the jet backgrounds can still be as large as 25% of the physical 77 background |]3"7fl , 
depending on the two-photon invariant mass. We adopted this figure as our estimate 
of the jet background that remains after all isolation cuts, applying it to all invariant 
masses. Finally, we assumed a 90% reconstruction efficiency for each photon as well as a 
further 7% loss, from isolation cuts, due to pileup, if the high luminosity option is used 



p7j . To obtain the number of background events, we used a mass bin of AM 77 = 3%M 77 . 
This is probably a conservative choice, given the energy resolution capability quoted in 
refs. [ p7| , p8[ . However, precise vertex identification is a very difficult task J7|, especially 



at high luminosity, and it is uncertain how accurately the vertex will be localized |37| . A 
conservative assumption on the vertex identification capability leads to an energy smearing 
such that our choice of AM 77 then corresponds, approximately, to a 2a bin width (see 



|37J). In the cases where the signal was broader than this mass bin, we summed the 
background over the width of the resonance. 

The signal event rate depends on all six parameters of the scalar sector as well as on 
the top quark mass and the choice of Model. The discussion of Sections 3 and 4 has shown 
that the most important parameter, besides M^o, is (3, by virtue of its effect on both the 
production rate and the two-photon decay width of the CP-odd Higgs. We will therefore 
present results fixing the other parameters at first. Nonetheless, we will discuss the effect 
of varying them and illustrate it by means of some examples. 

In Section 4 it was observed that the widths of the tt and scalar decay modes are 
generally large, so that the branching ratio to two photons is, most of the times, greatly 
reduced once these channels are open. Figures 7 and 8 reflect this behaviour. They display 
the production cross-section times the two-photon branching ratio of the A as a function 
of its mass for various values of (3. Fig. 7 corresponds to Set A of the parameters (see flilf)), 



while Fig. 8 corresponds to Set B (eq. (0))- We also plot, in each case, the same quantity 
for the Standard Higgs. In these calculations we included QCD corrections, as explained 
in Sections 3 and 4, but not cuts or efficiencies. The curves do not therefore represent the 
actual signal cross-sections, but, rather, they are shown for purposes of comparison. As 
expected from the discussion of the preceding sections, the rate increases with M A o until 
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the it or Zh threshold is met; beyond this, the rate falls due to the sharp decrease of the 77 
branching ratio. In Fig. 8, the cross-section rises again right before M A o = 2M t , reflecting 
the increased production rate at this mass. (The only exception to this general pattern is 
the case of large (3 in Model II; there, the quark loop diagrams responsible for both the A 
production from gluon fusion and the decay to two photons are dominated by the bottom, 
rather than the top, quark and so do not possess a maximum at M A o = 2M t .) Variation 
of M t , M H ±, Mjj, Mh may change the order in which the various thresholds appear, but 
has little effect otherwise. The dependence on the angle a is generally mild, except for its 
influence on the couplings of A to Zh and ZH, discussed at length in Section 4. Figure 9 
corresponds to Set C of parameters (see (p~5|) ) , an example where the suppression of the 
A°Zh coupling is particularly strong. 

Comparison with the Standard Higgs, in Figures 7-9, shows that in the so-called 
intermediate-mass range, the rates for the A are smaller, unless f3 ^ 40°. On the other 
hand, they can remain substantial for a much wider range, reflecting the fact that the 
WW and ZZ decay modes are effectively absent for the CP-odd scalar. 

In order to determine the statistical significance of a signal in this mode, we computed 
its rate using the cuts flTS] ) and the efficiencies quoted in the discussion of the backgrounds. 
We only accepted signals that consist of at least 40 events and defined the significance 
by the ratio S/y/B. In Figs. 10-12, we plot 5cr contours for integrated luminosities of 
10 ftr 1 and 100 fb _1 , and for each of Sets A, B and C. The contour corresponding to 
10 fb^ 1 can also be regarded as a 15a contour^] for an integrated luminosity of 100 fb _1 ; 
it thus serves also as an indication of the variation of the significance in the parameter 
space depicted. As explained in the Introduction, we restrict attention to the region 
of parameter space allowed by triviality constraints; this region is shown explicitly in 
Figs. 10-12. The significance plots show somewhat more quantitatively the behaviour 
expected on the basis of the results displayed in Figs. 7-9: For Sets A and B the opening 
of the it and Zh channels respectively, makes the detection of A virtually impossible 
beyond these thresholds, despite the small window around M A o = 2M t (Set B). On 
the other hand, the suppression of the strength of the A°Zh coupling leads to rather 
spectacular results for Set C. It is worth noting that this is not a case of fine-tuning: any 
a in the range —60° ^ a ^ —30° would have more or less a similar effect. 

Our comment on the dependence of the curves of Figs. 7-9 on the values of the other 
parameters of the model can help one visualize the variation of these significance contours 
with M t , M H ±, Mh, Mh and a: we expect the masses to play a role only insofar as they 

6 Here and in subsequent Sections we effectively scale the number of signal and background events with 
the integrated luminosity (except for the small pileup effect mentioned earlier). It is not absolutely clear, 
however, whether the assumed detector efficiency can be maintained in a high luminosity environment 

0. 
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set the various thresholds, and a to smoothly effect the transition from a situation like 
that of Fig. 11 to one like that of Fig. 12, as it gradually alters the factor of cos(/3 — a) 
governing the Zh branching ratio of the A . (The latter is relevant only if the Zh threshold 
comes before the it or W^H^ ones.) However, there is a different way in which all these 
parameters may affect our conclusions, and this is by changing the triviality bounds. The 
allowed area may be shifted to lower or higher (3, expanded or shrunk according to the 
precise values of the parameters chosen. This effect won't be too drastic if the scalars 
and top quark are not very light or very heavy; however, an example is shown in Fig. 13, 
where the region in which the A can be observed in this mode is rather narrow. We note 
the values chosen for this example for future reference: 

Set D : M H = 400 GeV, M h = 100 GeV, M H ± = 350 GeV, 

M t = 180 GeV, a = -75°. ^ ' 



In the MSSM, values of tan/3 > 1 are preferred |TJ]]. In this case, since a is not an 
independent parameter, cos(/3 — a) is constrained to be small, especially for large M40; 
thus the A°Zh coupling is suppressed, leading to a situation similar to that of Fig. 12b. 
However, the bias on tan (3 restricts us to a region of low production cross-section making 
this mode appear less valuable overall than in the Two-Doublet Model. 

As a conclusion, we may say that there is a substantial region of parameter space where 
the inclusive two-photon mode provides a clear signal for the A . This region extends from 
the lowest allowed (3 to (3 ~ 55 — 60°, with the signal fading as (3 increases, and from 
M^o ~ 40 GeV to the first threshold of a 'strong' decay (by which we mean it, W ± H T 
or Zh if not suppressed), with the signal improving with larger M^o. This is true for 
both Model I and Model II. From the examples considered it can be seen that Model II 
is slightly less favourable to high (3 values than Model I, but otherwise exhibits similar 
behaviour. These results, however, were based on rather optimistic assumptions about 
the ability of the detector to successfully reduce the jet backgrounds. This uncertainty 
has led many authors to seek alternative methods of detecting the scalar particles of the 
Standard Model and its extensions P,|9|,|3l| . We now turn to a method that has recently 
attracted a lot of attention. 



6 The /77X channel 

In the Standard Model, the Higgs boson can be produced in association with a W or a it 
pair through the following processes: 

qq -> W* -> WE (20) 
gg, qq -> ttH (21) 
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The leptonic decay of the W or the t/t gives rise to a final state that consists of two 
isolated photons from the Higgs and an isolated lepton (e or /j). Although the event rate 
is not high, the signal is rather clean, at least at low luminosity P,p|, ^0| , p| . [42 . 

For the CP-odd scalar of the Two-Doublet Model, (|20| ) gives a negligible rate, as 
stated earlier. Consequently, we concentrate on (f2~l~l). The final state from this process 



contains, besides the photons and the lepton, at least two jets. This can be used to remove 



most of the reducible backgrounds [40,42|, as will be explained below. 
The irreducible backgrounds arise from the processes 



99 
qq 



tt'yy 
ti'yy 



(22) 
(23) 



We did not perform an independent calculation of these backgrounds, but relied instead 
on the results of ref. [|12j for (^) []. The cuts assumed by these authors were the following: 



'•71 < 2.5 



1 77 

A J R( 7 i,7 2 ) > 0.4 



Pi 1 > 20 GeV, 
Ai2(Z,7) > 0.4. 



(24) 



Insisting that at least two jets be detectable in the final state results in a loss of rate that 
varies from 25% for M t = 100 GeV to 4% for M t = 160 GeV @. (In g| a jet was 
deemed detectable if it passed the same cuts as a photon or a lepton.) The background 
fl23|) was calculated in ref. and was found to amount to 50% of (^) at M 77 = 100 GeV 
climbing to 120% at M 77 = 180 GeV. We assumed these values in our estimate of the 
background. 

A potentially serious background could arise from the process gg —>■ 6677, where a 
B meson decays semileptonically. Since, however, the matrix element for this process 
is peaked when the b quarks are emitted along the beam direction, the cuts imposed 
on jets reduce this background to negligible levels |42|]. Other reducible backgrounds 
arise from events involving ti production in which one or two jets in the final state are 
misidentified as photons. In [T^] , the background from tig production was calculated and 
found negligible assuming a 7-jet rejection factor of 5 x 10~ 4 . A rough estimate of tt^ 
production was also reported and the resulting background declared to be small compared 
to the ttj'-f background, though not negligible. Another possible source of background is 
it production, where both photons are faked by jets from top decays. We assumed that 
these backgrounds can all be kept below the irreducible ^77 level. However, a detailed 
study is necessary in order to assess their exact magnitude at the LHC. 

7 These authors employed the HMRSB parton distribution functions p3| . We scaled their results by 
1.5 to account for the fact that EHLQ usually yield higher rates for gluon- initiated processes 0. 
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To obtain the background event number we used the reconstruction efficiency for 
photons and leptons (65% overall) quoted in |!7| and assumed, as before, a 3% resolution 
for the invariant mass of the photons. The background is a decreasing function of the 
top quark mass, due to the smaller parton luminosities at higher energies. Moreover, it 
depends indirectly on the mass of the charged Higgs: if the decay mode t — > H + b is 
kinematically allowed, then the leptonic (e and //) branching ratio of the top (or antitop) 
is reduced from the Standard Model value of 2/9 because H + decays to electrons or 
muons are very rare. This remark applies, of course, also to the signal cross-section, to 
the discussion of which we now turn. 

Fig. 14 shows the cross-section times branching ratio for the process pp — > it A —>■ 
/77A, as a function of M A o for various values of (3. The other parameters are chosen 
to belong to Set A (see (|i~3"D). The corresponding curve for the Standard Higgs is also 
shown. Cuts were not included in this calculation. The variation of the curves with 
the top quark mass is illustrated in Fig. 15, where M t = 180 GeV while the remaining 
Two-Doublet Model parameters are as in Set A. We shall refer to this choice as 'Set E'. 
We note that, due to the weak dependence of the production cross-section on the mass 
of the top, the rates are not too different for Sets A and E as long as M40 < 2M t . The 
production cross-section does not depend on parameters other than M A o, M t and /3, while 
the variation of the branching ratio with the other parameters has been discussed in 
Section 4. Consequently, we expect, as before, a sudden drop of the signal rate when the 
it or scalar thresholds appear, unless again there is a suppression of the A°Zh coupling. 
What is noteworthy in Figs. 14-15 is that for Higgs masses of less than about 160 GeV and 
for a large range of values of (3 the rate is well below its Standard Model value, especially 
in Model I. Considering the fact that a Standard Higgs of this mass can only give a few 
tens of events [f^ [4T|f41~| , we recognize that the prospects of detecting A in this mode are 
rather poor. In compensation, however, the mass range where this mode could be - even 
marginally - helpful may be much wider. 

In Figs. 16-17 we display significance contours for Sets A and E. Here the calculation 
of the signal was carried out including the cuts (|2"4[), reconstruction efficiencies and losses 
from insisting that at least two jets be detected in the final state. The fraction of signal 



lost due to this last requirement is the same as for the background (see ref. ^2[). The 



kinematical cuts result in a signal loss that varies from 75% at M40 = 30 GeV to about 
30% at M A o = 300 GeV (with a very soft dependence on the top quark mass). We 
required that a signal candidate should consist of more than 10 events. For low numbers 
of events, the formula for the significance used in Section 5 is not very accurate because 
it is based on Gaussian statistics. Since the production of 'events' is a Poisson process, 
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the probability that an apparent signal is in fact due to background fluctuations is [0] 



P = £ P n (B) (25) 

n=S+B 

where S, B are the numbers of expected signal and background events respectively and 
P n (B) is the probability that n background events occur, given by the Poisson distribution 

P n (B) = B n e- B /n\ . (26) 

(Since S + B is not an integer in general, an interpolation between the two integers closest 
to it is employed.) The probability P can then be converted into a significance. 

The contour for an integrated luminosity of 10 fb _1 is absent from Figs. 16-17; at such 
low luminosity the LHC cannot detect the CP-odd scalar of the two-doublet model in 
this mode at the 5cr confidence level. For Set E, Model I, even a luminosity of 100 fb _1 is 
inadequate. The main problem is not the background, which is small, especially for large 
M t , but rather the size of the signal. A reduced isolation radius of AR = 0.2 has the effect 
of boosting the signal by 20%. A similar answer was found for the background in ref. 
50]; an understanding of the effect of this change in AR on the reducible backgrounds is 



needed before deciding whether this procedure can enhance the significance of the signal. 

The chances of observing the A would be improved if the value of (3 were lower, since 
then the production would be raised. In Fig. 18 we show an example where low (3 values 
are allowed by triviality. (This can be achieved, roughly speaking, when \a\ is small, 
the top quark is light and the scalar masses moderate.) For this example, we chose the 
following values: 

Set F : Mh — 350 GeV, M h = 150 GeV, M H ± = 350 GeV, 

M t = 120 GeV, a = 0°. ^ ' 

We notice that for (3 ^ 30° a highly significant signal is possible, at least for Model II. 
Recall, though, our earlier comment that such low (3 values are in danger of being already 
ruled out by experiment. (Both the measurement of BR(b —>■ 57) and that of the Z —>hh 
decay width set lower bounds on (3; for charged Higgs masses of a few hundred GeV the 
bound from the latter is the strongest If this example represents the 'corner' of 

parameter space where the prospects of detection of the A in this mode are optimal, 
Set D (see fll9f) ) of parameter values is such that no point in the (M40 , 0) plane gives a 
5<t signal in either luminosity. 

Summarizing, we remark that the /77V mode may give a clear signal in only a narrow 
range of parameter space, due to the low A°tt production rate, which, for masses below 
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150 GeV approximately, is considerably smaller than the corresponding rate for the Stan- 
dard Higgs. Consequently, this mode was found to be useful only if j3 ^ 35°, a range that 
may easily be disallowed by triviality (or experiment). If the integrated luminosity at 
the LHC is 10 fb" 1 , then the A can be observed in this mode only in exceptional cases. 
Comparison with the inclusive two-photon mode reveals the latter, under the optimistic 
assumptions listed in Section 5, to be a more promising option. The Z77X mode can be 
a useful independent check, but enough events can only be accumulated if j3 is relatively 
small and the high luminosity option at the LHC is used. 



7 The Z77 channel 

In Section 4 we observed that the branching fraction of the A to the states containing 
other scalars (Zh, ZH, W ± H^) is generally large, and suggested that these modes may 
prove valuable in the search for the CP-odd Higgs, particularly for large M40 where the 
two-photon mode cannot be relied on. In this section we examine in detail the process 

A -> Zh -> / + Z~77 (28) 

where I stands for an electron or a muon. A similar discussion applies to the mode 
A — ► ZH. For defmiteness, we shall keep H heavy. Allowing for the possibility that H 
is light enough to have an appreciable branching ratio to two photons is straightforward 
to implement, but would not add to our understanding of the phenomenology of the A . 

The direct background to ( p8|) comes primarily from the process 

qq -> Z77 (29) 

It was found to be small, as we show below, so we did not consider necessary to perform 
the calculation of gluon fusion into Z77 (through the 'pentagon' diagram) which could in 
principle be comparable. There is another process that can interfere with (|29|), namely 
the standard production of a CP-even Higgs in association with a Z: 

qq^ Z* -> Zh (-> Z ll)- (30) 

This was also found to be very small, so we will neglect it in what follows. The calculation 
of (|29"D was done using spinor techniques |29|. The cross-section is singular when the 



photons are emitted along the beam axis or are soft. The following cuts were employed 
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to exclude such photons and simulate a realistic situation: 



|V' 7 | < 2.5 
Ai?(7i,7 2 ) >0.4 
AR(l,i) > 0.4 



Pi 1 > 20 GeV, 
AR(hJ 2 ) >0.4, 



(31) 



The isolation cut for photons and leptons has only a 5% effect on this calculation. We 
also included reconstruction efficiencies of 90% for each lepton and 85% for each photon 
after isolation cuts (see ref. pTfl). 

Fig. 19a shows a plot of the number of Z77 background events per GeV per LHC 
year (at a luminosity of 10 34 cm~ 2 s _1 ) as a function of the 77 invariant mass. For an 
invariant mass in the range 40 GeV ^ M 77 ^ 160 GeV (the interesting mass range for 
a Higgs h to decay to two photons), a 3% detector resolution will yield ~ 2 events of 
background. A similar conclusion is reached by looking at Fig. 19b, which plots events 
per GeV per LHC year versus the total final state invariant mass. For typical values of 
parameters that result in a signal of several events, the width of A is a few GeV. The 
background in such a bin is about 1 event. (The A can be considerably broader at high 
mass, but there the M Zll distribution dies out.) The h and A resonances will appear 
as bumps in the M 77 and M^ 77 distributions respectively, and it is the same events that 
will constitute both bumps. This is not the case with the background, however. The 
background events with a two-photon invariant mass M 77 within 3% of Mh will have a 
broad distribution in M^ 77 and will not, in general, be concentrated in the vicinity of the 
A resonance. This observation enables us to effectively eliminate the background. For 
example, Higgs masses of M h = 100 GeV, M A o = 300 GeV may typically give a signal of 
a few tens of events. Of the ~ 2 background events that have M 77 inside a 3% window 
around 100 GeV, less than 0.1 lie within 5 GeV of M^o. Hence it becomes plain that the 
background that simultaneously emulates both h and A is truly negligible. 

QCD corrections to process ( p9|) have not been calculated. We also did not estimate 
the size of reducible backgrounds from jet events (Zjj, Zjj). It is assumed that, with the 
efficiencies and rejection capabilities mentioned in Section 5 [571, these backgrounds will 
be brought to below the level of the direct Z77 production. But in view of the arguments 
of the preceding paragraph, even if the reducible backgrounds were to dominate by as 
much as an order of magnitude over the irreducible Z77 production, they would still pose 
no significant problems. 

It is evident, therefore, that the observability of the A in this mode depends only 
on whether a sufficient number of signal events can be accumulated. The magnitude of 
the light neutral Higgs branching ratio to two photons, BR(h — > 77), is thus of crucial 
importance. This branching ratio has a complicated dependence on the various masses 
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and angles. We will indicate very roughly what this dependence is so as to make the 
pattern of our results more intelligible. 

A Higgs h in the intermediate- mass range 40 GeV ^ ^160 GeV, decays mostly 
to a bb pair. Therefore, to a good approximation, 

From the hbb couplings displayed in Table 1, it follows that 

r (^^)~{ cos 2 2a / sin ^ ( ^Vl ^ 

v ; \ sin 2 a/ cos 2 /3 (Model II) v ; 

The partial width T(h — > 77) is more complicated. The light neutral Higgs h decays 
to two photons through a loop in which a W, an H + or a charged fermion flows. For 
moderate values of (3 the contribution of the W is the largest, although, in some cases, 
the top quark loop is of comparable magnitude. The hWW coupling (see Table 1) is 
proportional to sin(/3 — a) so we may write, very roughly, 

nn/i \ f sin 2 (/5 - a) sin 2 /5/cos 2 a (Model I) , nA . 

BR{k ~> 77) ~ 1 sin 2 (/3 - a) cos 2 (3/ sin 2 a (Model II) (34) 

provided none of the sines and cosines that appear in these formulae are too small. An- 
other feature of the 77 branching ratio of the scalar h is that it becomes largest for « 
120-130 GeV and has a very soft dependence on M H ± and M t . 

Our prime motivation for considering this decay mode was to explore regions of pa- 
rameter space where the two-photon channel gives a signal too faint to observe. In Fig. 20 
we show Z77 event contours for Set D of parameters (see (|i~9|)). More than 10 events are 
expected for points inside these contours. The solid curve (corresponding to an integrated 
luminosity of L = 10 fb _1 ) can also be interpreted as a 100-event contour for L = 100 ftT 1 . 
The branching ratio for the Z decay into electrons or muons, as well as the kinematical 
cuts (|5ID and the lepton and photon reconstruction efficiencies, were included in the cal- 
culation of the signal. QCD corrections to the decays of h were taken into account in a 
way similar to that described in Section 4 for the A 0. The large value of \a\ (see (|19|)) 
results in a huge rate in the case of Model I, where, at L — 100 fb^ 1 , a clear signal is 
obtained for almost the entire parameter space depicted. For Model II the results are 
more modest, but comparison with Fig. 13 shows that, even in this case, the ^77 signal 



3 In addition, the known corrections to the two-gluon mode ]34f] were also included. 
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is observable in a region where the two-photon mode can hardly be helpful. Choosing \a\ 
to be small would instead produce a large rate in Model II. This is illustrated in Fig. 21. 
The only difference relative to Set D is the value of a which in this example is a = 0°. 
We shall refer to this choice as 'Set G'. Here too, the inclusive two-photon mode fails to 
give an observable signal for A masses beyond the threshold for the Zh decay. (We omit 
plots for this, however.) 

In Section 5 we remarked that the 77 signal is weak for large M A o , when other channels 
open up, and for large /3, when the production rate (as well as the two-photon branching 
ratio in Model II) is small. From Q53|) it follows, however, that, depending on a and the 
choice of Model, the light CP-even Higgs branching ratio to two photons may be enhanced 
at large /3 and perhaps even compensate for the loss in production rate. An example where 
this happens is given in Fig. 22. The parameters belong to Set C, which represents the 
case where the A°Zh coupling is weak enough to allow a sizeable 77 branching ratio to 
persist up to the top threshold (see Fig. 12). Indeed, the Zh mode can't give a good 
signal beyond this threshold either, but, at high luminosity, it covers some of the large-/? 
regions undetectable by the two-photon mode, while being able to confirm the latter for 
smaller (3 and M z + M h < M A o < 2M t . 

In order to convey an idea on the dependence of our results on a, we present in Fig. 23 
similar 10-event contours in the (a,/3) space, now keeping M A o fixed: 

Set H : M H = 400 GeV, M h = 100 GeV, M H ± = 350 GeV, 

M A o = 280 GeV, M t =150 GeV. ^ ' 

The precise value of M40 has no effect on these curves, as long as it is below the W ± H T or 
tt threshold. The different behaviour of the two Models for large and small \a\ is evident 
from the figure. We observe that large numbers of events are expected in a considerable 
region of parameter space. The only area left unexplored corresponds to a > and 
large (3. 

The situation is not as rosy if M40 > 2M t or M40 > M w + M H ± , in which case the 
branching ratio to Zh is smaller. A clear signal of a few tens of events is then guaranteed 
in Model I (resp. Model II) only for large (resp. small) |a|. 

In all figures shown the value Mh = 100 GeV was assumed. We remarked earlier 
that the branching ratio of h to two photons increases with Mh until approximately 
Mh ~ 120 GeV and then falls again. The corresponding Z77 event rates would of course 
follow the same trend. 

In the MSSM the decay sequence A — > Zh — > / + /~77 discussed in this section has 
not been considered, presumably because it gives too low a rate. We have mentioned, in 
Section 5, that the bias on tan/3 tends to suppress the A°Zh coupling in the kinematical 
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region where the decay A — > Z/i is possible. Consequently, the branching ratio of this 
decay is only sizeable when tan (3 is not very large and Mz + < M40 < 2M t . How- 
ever, in this range, the branching ratio of h to two photons does not exceed 10~ 3 (except 
for large M t ) [111. In contrast, in the context of the non-supersymmetric Two-Doublet 



Model, it can reach the level of a few percent (and even higher in exceptional cases). 
The different behaviour can be traced to the fact that a is an independent parameter in 
the Two-Doublet Model and may assume values such as to enhance this branching ratio. 
The decay A — > Zh in the MSSM has received some consideration in conjunction with 



the subsequent r + r [44] or bb uM decays of the light Higgs h, with some encouraging 



results in the former case. The r + r~ branching ratio of an intermediate-mass h in the 
Two-Doublet Model, as well as in the MSSM, is about 10%, almost independently of a 
and (3. Consequently we may expect this mode to be helpful also in the Two-Doublet 
Model; in particular, it may enable us to explore the region for a > where the Z77 
signal is weak. A detailed study of the signal and backgrounds at the LHC is needed in 
order to decide this issue. 

To summarize, we have found that the process A — ■> Zh — > provides an 

excellent way of simultaneously detecting the A and the light CP-even scalar h, provided 
the latter is in the intermediate mass-region 40 GeV ^ Mh ^ 160 GeV. This is true even 
at a yearly luminosity of 10 fb -1 at the LHC. Tens or hundreds of signal events can be 
obtained over a negligible background in a fairly large region of parameter space, especially 
when I a I is large (in Model I) or small (in Model II). This mode only fares poorly when 
a is positive and moderate and (3 large, or when the ti or W^H^ channels are available 
for the decay of A . If the latter is true, detection is harder but not impossible. 

The utility of the Zh channel extends also to larger masses for the CP-even Higgs h. If 
Mh ^ 160 GeV, then h can be detected in its 'gold-plated' decay mode h — > ZZ or ZZ* — > 
l + l~l + l~, with I = e,/i. We have estimated that the process A — > Zh followed by the 
above decay sequence may yield several events provided M^o <2M t . The last requirement 
shows that the kinematical region where this six-lepton signal can be observed is rather 
narrow, unless the top quark is relatively heavy (M t ^ 150 GeV). Nevertheless, this may 
be the only way of detecting the CP-odd scalar in this region of parameter space. 



8 Conclusions 

In this paper we examined the phenomenology of the CP-odd scalar A of a two-Higgs- 
doublet model with an exact discrete symmetry. We explored the parameter space defined 
by the triviality bounds of ref. || in order to identify the regions where the A can be 
detected at the LHC. We studied three different signals of the A . The inclusive two- 
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photon mode can provide an observable signal if j3 ^ 55-60° and provided the decays of 
the A to it, W ± H Zf , Zh are kinematically forbidden or (in the case of Zh) suppressed. 
The lower end of the A mass range where detection in this mode is possible depends on 
(3 and is raised as (3 increases. Models with small (3 lead to larger signals because the 
production rate of A is effectively proportional to cot 2 (3, at least for moderate values of 
(3. This is also true for the associated it A production, which was found capable of giving 
a detectable signal in the l^'yyX mode if (3 ^ 35-40° (and for roughly the same range of 
M^o as in the inclusive two-photon case). However, this range is mostly disallowed by 
triviality bounds (as well as disfavoured by current experimental data). If the integrated 
luminosity at the LHC is 10 fb _1 , this mode can only be helpful in exceptional cases. In 
contrast, the process A — > Zh, followed by the two-photon decay of the light CP-even 
Higgs h, can provide a very clear signal for a substantial region of parameter space, as 
long as h is in the so-called intermediate- mass range (40 GeV ^ Mh ^ 160 GeV). (For 
heavier h, its decay to four leptons through a Z pair may give rise to several events if 
the masses are such that Mz + < M^o < 2M t .) This mode is largely complementary 
to the 77 mode since it covers regions of high M^o as well as (in some instances) high (3. 
Under certain circumstances, it may also confirm observation of the A in the two-photon 
channel. Hundreds of events can be accumulated in this mode if |a| is large in Model I 
or small in Model II. 

If M^o > 2M t or M A o > Mw + Mff± , then the branching ratio to Zh is reduced and 
our ability to detect the CP-odd scalar in this mode is compromised. Although for some 
values of a and (3 this mode would still be a prime option for the direct observation of 
A , it may be worthwhile investigating other processes, such as the it A production with 
the subsequent decay of A to it |R|. Another region that, according to our conclusions, 



escapes observation at the LHC, is the low-M^o, large-/? region (essentially regardless of 
the values of the other parameters). It may be possible to explore part of this region at 
LEP-II, particularly since the disadvantage of insufficient production at large (3 will not 
be present there. 

It would be interesting to perform similar studies for the other scalars of the Two- 
Doublet Model, in order to determine how much of the parameter space of such models 
can be explored at future colliders. Work in this direction is under progress. 
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Figure captions 



1. Production cross-section of the CP-odd Higgs A through gluon fusion as a function 
of its mass for (a) Model I and (b) Model II. Results are shown for f3 = 15°, 45°, 70° 
and 88° and for M t = 150 GeV. 

2. Production cross-section of the A through bb fusion as a function of its mass for 
(a) Model I and (b) Model II and for the same set of values of (3 as in Fig. I. 

3. Cross-section for the associated it A production as a function of the A mass for (a) 
M t = 120 GeV and (b) M t = 180 GeV and for the same set of values of (3 as in Fig. 1. 
There is no distinction between Models I and II in this case. The corresponding 
cross-section for the Standard Higgs is also shown for comparison. 

4. Total production cross-section of the A as a function of its mass for (a) Model I 
and (b) Model II and for the same set of values of (3 as in Fig. 1. The corresponding 
cross-section for the Standard Higgs (not including gauge boson fusion production) 
is also shown. The top quark mass is taken to be 150 GeV. 

5. Branching ratios of the A as a function of its mass for the following choices of pa- 
rameters: (a) M H = 400 GeV, M h = 260 GeV, M#± = 350 GeV, M t = 120GeV, a = 
30°, (3 = 30°, Model I; (b) same for Model II; (c) same as (a) except f3 = 60°; (d) 
same as (b) except (3 = 60°. 

6. Branching ratios of the A as a function of its mass for the following choices of pa- 
rameters: (a) M H = 400 GeV, M h = 100GeV, M H ±=350 GeV, M t = 180GeV, a = 
30°, f3 = A0°, Model I; (b) same for Model II; (c) same as (a) except a = -45°; (d) 
same as (b) except a = —45°. 

7. Production cross-section times branching ratio to two photons for the A as a func- 
tion of its mass for (a) Model I and (b) Model II and for the following set of param- 
eters: M H = 400 GeV, M h = 260 GeV, M H ± = 350 GeV, M t = 120 GeV, a = 30°. 
Curves are shown for (3 = 15°, 45°, 70° and 88° as well as for the Standard Higgs. 

8. Same as Fig. 7 but for the following values of parameters: Mh = 400 GeV, Mh = 
100 GeV, M^±=350 GeV, M t = 180 GeV, a = 30°. 

9. Same as Fig. 7 but for the following values of parameters: M H = 400 GeV, M h = 
100 GeV, M H ±=350 GeV, M t = 180GeV, a = -45°. 

10. Significance contours for the inclusive two-photon signal in the region of the (M^o , f3) 
plane allowed by triviality for (a) Model I and (b) Model II. The other parameters 
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are as in Fig. 7. Contours are shown for integrated luminosities of 10 fb and 
100 fb -1 . The statistical significance is greater than 5cr in the interior of these 
contours. 

11. Same as Fig. 10 but for the set of parameters of Fig. 8. 

12. Same as Fig. 10 but for the set of parameters of Fig. 9. 

13. Same as Fig. 10 but for the following choice of parameters: M# = 400 GeV, Mh = 
100 GeV, M H ±=350 GeV, Mt = 180GeV, a = -75°. 

14. Associated it A production cross-section times branching ratio to an l ±/ ~fj final state 
as a function of M^o. This branching ratio is the product of the A — > 77 branching 
fraction and the probability that at least one of the t and t decays leptonically. 
Results are shown for (a) Model I and (b) Model II and for (3 = 15°, 45°, 70° and 88°. 
The remaining parameters are as follows: M H = 400 GeV, M h = 260 GeV, M H ± = 
350 GeV, M t = 120 GeV, a = 30°. The corresponding quantity for the Standard 
Higgs is also plotted for comparison. 

15. Same as Fig. 14 except now M t = 180 GeV. 

16. Significance contours for the /77X mode in the region of the (M A o , (3) plane allowed 
by triviality for (a) Model I and (b) Model II. The other parameters are as in Fig. 14. 
No point in the region depicted can give a 5a signal if a yearly luminosity of 10 fb -1 
is employed. 

17. Same as Fig. 16, but for the set of parameters of Fig. 15. 

18. Same as Fig. 16, but for the following values of parameters: M# = 350 GeV, M^ = 
150 GeV, M H ±=350 GeV, M i = 120GeV, a = 0°. 

19. Continuum Z77 background events per GeV per 100 fb -1 as a function of (a) the 
two-photon invariant mass M 77 and (b) the total final state invariant mass M Zll . 

20. 10-event contours for the Z77 signal in the region of the (M A o , (3) plane allowed by 
triviality for (a) Model I and (b) Model II. The slice of parameter space shown corre- 
sponds to the following set of parameters: M H = A00 GeV, M h = 100 GeV, M H ± = 
350 GeV, M t = 180 GeV, a = —75°. Contours are shown for luminosities of 10 fb -1 
and 100 fb -1 . More than 10 events are expected in the interior of these contours. 

21. Same as Fig. 20 except now a = 0°. 

22. Same as Fig. 20 except now a = —45°. 
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23. 10-event contours for the Z77 signal in the region of the (a , /3) plane allowed 
by triviality for (a) Model I and (b) Model II. The remaining parameters take 
the following values: M H = 400 GeV, M h = 100 GeV, M H ± = 350 GeV, M t = 
150 GeV, M A o = 280 GeV. 
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